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Abstract 



We consider the lifetime of quasi-particles in a d-wave superconductor due 
to scattering from antiferromagnetic spin- fluctuations, and explicitly separate 
the contribution from Umklapp processes which determines the electrical con- 
ductivity. Results for the temperature dependence of the total scattering rate 
and the Umklapp scattering rate are compared with relaxation rates obtained 
from thermal and microwave conductivity measurements, respectively. 
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It is now widely accepted that the superconducting state of the high T c cuprates is 
characterized by a rf a .2_ y 2-order parameter Furthermore, various transport properties 

such as NMR relaxation times and the microwave conductivity H have been calculated 
within a quasi-particle BCS-RPA framework which phenomenologically takes into account 
the interplay of both the d x 2_ y 2 pairing and the antiferromagnetic spin correlations. In 
particular, this approach has previously provided reasonable fits of the temperature de- 
pendence of the NMR T\ and T 2 relaxation times in optimally doped YBa 2 Cu307_,5 H0. 
With recent advances in the preparation of YBCO crystals HH, microwave measurements 
|T0| at a variety of higher frequencies, and K xy thermal conductivity measurements JTTJ , one 
has new information on the quasi-particle lifetime involved in the electrical conductivity and 
the mean-free path associated with the thermal conductivity. Several puzzles have presented 
themselves. First, microwave measurements appear to indicate that the relaxation rate, long 
known to collapse rapidly below the critical temperature f[2-[14}|, varies as T 4 , or possibly 



exponentially, at temperatures T -C T c . On the other hand, the theory of the quasi-particle 
relaxation rate in a c?-wave superconductor due to electron-electron interactions predicts a 
T 3 variation above a scale set by impurity scattering J5|, [T3| , P~5| ] . A solution to this discrepancy 



has recently been suggested by Walker and Smith jT7|, who noted that in the absence of 
disorder, the transport lifetime entering the microwave conductivity arises solely from Umk- 
lapp processes, which are always gapped below T c even in a d-wave superconductor, leading 
to an exponential decay 1/tjj ~ exp(— Ajj/T). Secondly, while early thermal conductivity 
experiments extracted quasi-particle mean-free paths in apparent agreement with microwave 
measurements, recent measurements indicate that on cleaner crystals, where the mean-free 
path due to inelastic scattering can be probed to lower temperatures, the extracted mean- 
free path is shorter than in the microwave case ||11|| . Furthermore, the inelastic contribution 
to the inverse thermal mean free path appears to follow a T 3 temperature dependence at 
low temperatures (see below). 



Motivated by the proposal of Walker and Smith [17j regarding the role of Umklapp pro- 
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cesses, we investigate whether the same type of spin-fluctuation interaction used to fit the 
NMR data can also provide a quantitative fit to the lifetimes extracted from the microwave 
and thermal conductivities. In particular, we wish to understand if the discrepancies dis- 
cussed above can be resolved by recognizing that normal electron-electron interactions con- 
serve the total momentum and therefore cannot degrade the electric current, whereas they 
can degrade the heat current We therefore calculate, within a BCS model in which the 
spin-fluctuation interaction is accounted for in the random phase approximation (RPA), the 
total quasi-particle lifetime and the Umklapp-only quasi-particle lifetime and compare these 
with the K xy thermal conductivity and microwave conductivity measurements, respectively. 
We find that the inelastic lifetimes extracted from thermal transport and charge transport 
can be accounted for within a BCS-RPA quasi-particle framework. Furthermore, the param- 
eter values which are required are consistent with what is known about the Fermi surface 
P"9^PH and spin- fluctuations from other measurements ||[7j . The analysis is based, however, 
on the subtraction of a constant relaxation rate from both quantities, an approach which 
must be justified within microscopic theory; we present a brief discussion of these issues 
and the light shed by our analysis on the behavior of the low-energy impurity scattering 
amplitude. 

In the following, we will assume that we have BCS quasi-particles with a d x 2_ y 2 gap 



moving on a 2D square lattice with near- and next-near-neighbor one-electron hopping 
matrix elements t and t', respectively. The effective electron-electron interaction will be 
modeled by the spin-fluctuation form 



Here, U and the chemical potential, fi, which sets the average site occupation, (n), control 
the strength of the antiferromagnetic spin- fluctuations, and Xo(q> u ) is the BCS susceptibility 
given by 



A k (T) = -A (T)(cosk x - cosky) 



(1) 
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where E p = ^jep + A p and e p = — 2t(cosp a; +cosp 2/ )— 4i'cosp a; cosp 2/ — fi. The band parameters 
are taken to be if ft = —0.35 with a site filling of (n^ + n^i) = 0.85, and the temperature 
dependence of the gap has been parametrized as 



A (T) = ^tanh( av /r c /r-l). 



(4) 



We have used a = 3.0 to reproduce the rapid onset of the gap below T c which is observed 
experimentally, and set 2A /fc B T c = 5. Throughout this paper, we use t = 1 as our energy 
scale and take U/t = 2.2; this is similar to the U = 2.0 used in previous calculations |J of 
the NMR T\ and T2 relaxation rates, where an RPA form for x(q, iS) similar to Eq. (|^) was 
employed, but a simple (if = 0) tight-binding band assumed. 

Using the effective interaction, Eq. (|2|), the Fermi Golden Rule expression for the quasi- 



particle lifetime is [15 
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where n(u) and /(u;) are the usual Bose and Fermi factors. Here the first term arises from 
quasi-particle scattering resulting from the absorption and emission of a spin-fluctuation. 
The second term corresponds to a process in which a quasi-particle recombines with another 
quasi-particle to form a pair and the excess energy is emitted as a spin-fluctuation. 

An example of the type of process which this includes is schematically illustrated in 
Fig. 1(a). Here, a quasi-particle with momentum p scatters off of another quasi-particle 
with momentum k leading to an intermediate state containing quasi-particles of momentum 
p' and k'. Momentum conservation on the lattice implies that 



p + k = p' + k' + G (6) 

where G is a reciprocal lattice vector and all other momenta are confined to the first Brillouin 
zone. Since E^/ = E p _ p r +k _ G = -Ep_p' + k, we have the expression for xoi.*!^) given by 
Eq. (|3|) with q = p — p' and with the k sum extending over the entire first Brillouin zone. 
If, however, we had not allowed any reciprocal lattice vectors except G = (0, 0), the k sum 
would, for given values of p and p', extend only over that part of the first Brillouin zone for 
which k' = p — p' + k remains inside the first Brillouin zone. This restriction of G = (0, 0) 
eliminates all of the Umklapp scattering processes, allowing only "normal" processes. In this 
way, one can introduce a "normal" scattering lifetime T^ 1 (p,cu). This is obtained by using 
an interaction with the same denominator as in Eq. (Q), but with Xo(q, to) in the numerator 
replaced by the imaginary part of Eq. (|3]) with q = p — p' for which the sum over k is 
restricted such that p — p' + k remains in the first Brillouin zone. With this restriction, the 
resulting "normal", non- Umklapp, part of Xo(q, to) depends separately on p and p' rather 
than only on the difference of p — p'. Having calculated both r _1 (p, ui) and r^ 1 (p,u), one 
can obtain a lifetime which contains only the Umklapp processes by subtraction 

^(p^) = T~ x {p,u) - rj^fau). (7) 

At temperatures which are small compared with the maximum gap, the thermally ex- 
cited quasi-particles occupy states close to the d x 2_ y 2 gap nodes on the Fermi surface. When 
two of these thermal quasi-particles scatter, the outgoing states must also occupy states 
near the nodal region in order to conserve energy. For the Fermi surface shown in Fig. 1(b), 
these scattering processes must have G = (0, 0) and only contribute to the "normal" pro- 
cesses. Umklapp scattering processes are also possible, but these processes are exponentially 
suppressed at low temperatures. Following Walker and Smith [T7|], the reason for this is il- 
lustrated in Fig. 1(b). Here, we have set the initial quasi-particle momentum p to a node 
and numerically searched for the lowest energy quasi-particle state k such that momentum is 
conserved as in Eq. (||) with the constraint that G = (0, 0). For t'/t = —0.35 and (n) = 0.85 
corresponding to the Fermi surface shown, the momentum k is the lowest energy state for 
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which a quasi-particle can have an Umklapp scattering with the nodal quasi-particle p. Since 
k has moved away from the nodal region, A*, is now finite. When k^T is less than A^, the 
probability of finding such a quasi-particle varies as e _Ak//fcsT and the Umklapp scattering is 
suppressed for temperatures below Ak. We will call this special value the Umklapp gap, A^. 
It depends upon the shape of the Fermi surface which is set by t' j t and (n) . Figure 2 shows 
the Umklapp g function of the strength of the t'/t for a fixed filling of (n) = 0.85. It 

is easily seen that an Umklapp gap exists even when t'/t = and has a magnitude of about 
0.07t for the band with t'/t = —0.35. We note further that this corresponds to a "Dirac 
cone anisotropy ratio" Vp/v% ~ 12, which is consistent with the value determined by low 



temperature thermal conductivity measurements [|TI| . 

We now compare the absolute values of the relaxation times in question. In Fig. 3(a), 
we have plotted h/k B T c Tjj and 7i//cbT c t versus T/T c for the parameters discussed above. In 
order to compare recent thermal conductivity measurements with the total scattering rate, 
we need to specify a Fermi velocity. For the t-t' dispersion relation one has vp — 2.16ta 



and taking T c = O.lt = 90K, this gives Vp ~ 10 7 cm/sec. Zhang et. al. [IT] have used 
an effective Fermi velocity rjVf ~ 10 7 cm/sec with rj = 0.6 and vp = 1.7 x 10 7 cm/sec. 
Alternatively, Chiao et. al. [0| have suggested that vp = 2.5 x 10 7 cm/sec. In the fit of 



h/kBT c r(T) shown as the solid line in Fig. 4(a), we have used vp = 3x 10 7 cm/sec to convert 
the thermal mean-free path data of Ref. ]TTJ] to a quasi-particle relaxation rate shown as 
the open symbols. Beyond the choice of Vp, the determination of the quasi-particle lifetime 
from the thermal conductivity measurements is complicated by the phonon contribution to 
the heat current, and necessitates certain assumptions regarding the nature of quasi-particle 
scattering in the vortex state which are still controversial. However, we note that the 



magnitude of the theoretically determined total scattering rate of h/r(T c ) ~ 3ksT c at the 



node is in rough agreement with that determined by ARPES f22" |, while K/tu(T c ) ~ ksT c 
is consistent with optical data. In Fig 3(b), we show the low temperature range of the 
theoretical predictions, along with the Hosseini et. al. JTTJ 1/r extracted from microwave 



data using a two-fluid analysis, shown as open diamonds. Here one sees that taking U/t = 2 



provides a reasonable fit to the magnitude of the low temperature microwave scattering 
lifetime and the exponential decay of the Umklapp processes also appears to provide a 
sensible explanation for the observed temperature dependence at low temperatures. 

Clearly, there are a number of parameters in this calculation, and the results do vary 
with their values. Here, we have selected values for t'/t and the filling (n) to give a Fermi 
surface which resembles that seen in ARPES data and various band structure calculations 
20fl for the d x 2_ y 2 — pa band of YBCO. Specifically, the Fermi surface, as seen in Fig. 1(b), 
has a diagonal nodal k-spacing of (0.8, 0.8)7r/2a, and fixes the size of the Umklapp gap, Ajj, 
which determines how r^ 1 (p, o>) decreases below the total scattering lifetime. Note that for 
the Fermi surface pictured in Fig. 1(b), 2A / ksT c = 5 corresponds to a maximum of the 
gap on the Fermi surface near (ir, 0) of max[2Ao(p)//c_BT c ] = 4.7. The choice of U/t = 2.2 
was made consistent with Refs. [[|[|| and because it provided a reasonable fit to the absolute 
magnitude of the microwave scattering lifetime, as discussed below. 

It is important to ask the extent to which the data support the fundamental prediction of 
this work, that the relaxation rate entering the thermal conductivity data should be the total 
quasi-particle relaxation rate, which varies as (T/T c ) 3 due to spin-fluctuation scattering and 
recombination processes, in contrast to the activated Umklapp scattering rate which should 
enter the microwave conductivity. Is Fig. 3 sufficient to conclude that these assertions are 
correct? In order to examine the true low-T asymptotic forms, we have displayed in Fig. 4 
the same curves and data as in Fig. 3 on a log-log plot. ^From the insert to the figure, it is 
clear that the low-T behavior of the two experimental relaxation rates is quite different, but 
although the thermal conductivity data is roughly consistent with the T 3 result expected 
from the spin-fluctuation theory, the microwave data appear to lie considerably above the 
activated prediction at the lowest temperatures, a behavior hidden by the linear plot in Fig. 
3. In the main part of the figure, we show the raw experimental data with no constant 
subtracted, together with the inelastic rates added to a constant l/r determined by the 
experimental data at low T. 

The most natural interpretation of the residual scattering rates 1/tq is in terms of dis- 



order. It is striking that the 1/tq extracted from the thermal conductivity measurement 
appears to be significantly larger than that extracted from the microwave measurements, 
since both measurements were performed on nominally equivalent BaZrOs crucible-grown 
samples. The most popular treatment of disorder in the cuprates is the self-consistent 
t-matrix approximation (SCTMA) describing strong scattering by in-plane near-unitarity 
limit defects. This theory predicts a strong energy dependence of the low-temperature elas- 
tic scattering rate, however, inconsistent with the roughly constant behavior observed here 
15fl . This difficulty was remedied to some extent in Ref . |23| by considering additional scat- 



tering from spatial fluctuations of the order parameter. Even in this approach, however, bare 
scatterers are considered to be isotropic, leading to a vanishing of impurity vertex corrections 
in local quantities; this implies the same relaxation time for both thermal and microwave 
conductivities. If the isotropic scattering assumption is relaxed, however, as considered by 



Durst and Lee [[OJ, the uo — > microwave conductivity is found to be strongly enhanced 
relative to the thermal conductivity, which is unrenormalized. This may be interpreted as a 
reduction of the scattering rate in the microwave case, but a full frequency-dependent cal- 
culation and proof that a Lorentzian conductivity obtains in analogy to Ref. |J is required 
before one may conclude that the observed discrepancy is caused by anisotropic impurity 
scattering. 

To summarize, using a d-wave BCS framework, we have taken a simple spin-fluctuation 
interaction, which has been used to describe the temperature dependence of the NMR T\ 
and T 2 spin relaxation times in YBCO below T c , and have shown that it can also describe the 
temperature dependence of both the inelastic charge and heat transport times below T c . The 
contribution to the quasi-particle lifetime of normal spin-fluctuation scattering processes, 
which can relax heat currents but not charge currents, as well as Umklapp processes, have 
been calculated separately. The difference in the two contributions was then shown to 
account for the difference in experimentally extracted lifetimes in recent microwave and heat 
conductivity measurements on ultraclean samples of YBCO. Both the predicted temperature 
dependences in a spin-fluctuation mediated interaction model, the Umklapp scattering rates 
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1/tjj ~ exp —Au/T and the total scattering rate 1/r ~ T 3 , appear to agree with the data, 
as do the predicted absolute magnitudes of the two rates. The exponential form of 1/tjj 
implies a sensitivity to the Fermi surface shape, which is expected to lead to a significant 
increase in 1/tu at low temperatures upon underdoping. While our results depend upon 
the band structure and interaction parameters, we have selected these to be consistent with 
Fermi surface data and NMR measurements, and conclude that the results for the microwave 
and thermal conductivity inelastic quasi-particle lifetimes are consistent with an effective 
electron-electron interaction which is mediated by spin-fluctuations. 
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FIGURES 
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FIG. 1. (a) A schematic representation of the scattering process, (b) The non- interacting 
Fermi surface with t'/t = —0.35 at a filling of (n) = 0.85. The vectors shown correspond to an 
Umklapp process in which p + k = p 1 + k' + (0, 2ir/a). 
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FIG. 2. The size of the Umklapp gap Ajj/t (solid curve) and the vpjv2 ratio (dashed curve) 
versus the ratio of the next-near-neighbor hopping t' to the near-neighbor hopping t. 
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FIG. 3. (a) The normalized YBCO quasi-particle scattering rates 7i/(t/cbT c ) and h/ '(t[/&bT c ) 
versus reduced temperature T/T c . The experimental points for K/ (rksTc) were obtained from the 
thermal conductivity mean- free path measurements by Zhang, et. al. [|ll] using up = 3x 10 7 cm/sec 
and are shown as the open circles and squares. The solid line is the total quasi-particle lifetime 
obtained from Eq. (§) on a 512 x 512 lattice with U/t = 2.2, T c = O.li, t'/t = -0.35, 2A /k B T c = 5, 
and (n) = 0.85. (b) Same as in (a) for reduced temperatures less than 0.5 with the data of Hosseini, 
et. al. [^] included as open diamonds. The solid line corresponds to the total lifetime calculated 
as in (a), while the dashed line represents the lifetime due to Umklapp processes only, Eq. ([?]). 
Note a constant term has been subtracted from the experimental data. 
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FIG. 4. Log-log plot of theoretical predictions for normalized total and Umklapp scattering 
rates 7i/(rfesT c ) and h/(rukBT c ) obtained from Eq. (||) (solid lines), same parameters as Fig. ||. 
Data from Refs. ||Tl| and |fTcf| as in Fig. §. Dashed curves: same theoretical predictions for 
h/(Tk B T c ) and fi/{Tjjk B T c ) with constants l/r = 5.0 x 10 10 s _1 and 3.3 xlO 11 s _1 added to 1/tu 
and 1/t, respectively. Insert: same as main panel but with 1/tq subtracted from data. 
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